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Abstract: We construct a set of parton distribution functions (PDFs) in which fixed-order 
NLO and NNLO calculations are supplemented with soft-gluon (threshold) resummation 
up to NLL and NNLL accuracy respectively, suitable for use in conjunction with any QCD 
calculation in which threshold resummation is included at the level of partonic cross sections. 
These resummed PDF sets, based on the NNPDF3.0 analysis, are extracted from deep- 
inelastic scattering, Drell-Yan, and top quark pair production data, for which resummed 
calculations can be consistently used. We find that, close to threshold, the inclusion of 
resummed PDFs can partially compensate the enhancement in resummed matrix elements, 
leading to resummed hadronic cross-sections closer to the fixed-order calculations. On the 
other hand, far from threshold, resummed PDFs reduce to their fixed-order counterparts. 
Our results demonstrate the need for a consistent use of resummed PDFs in resummed 
calculations. 
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1 Introduction 

The accurate determination of the parton distribution functions (PDFs) of the proton is 
an essential ingredient of the LHC physics program [1-5]. In order to reduce theoretical 
uncertainties, it is crucial to incorporate in global PDF fits higher-order perturbative QCD 
corrections, both to the hard partonic cross sections and to the parton evolution. While 
recent progress in fixed-order NLO (see e.g. [6] for a recent review), NNLO (e.g. [7-23]) 
and even N^LO [24] calculations for different processes in hadron-hadron collisions has 
been impressive, it is also well-known that fixed-order perturbative calculations display 
classes of logarithmic contributions that become large in some kinematic regions, thus 
spoiling the perturbative expansion in the strong coupling constant The importance 
of these contributions varies significantly with both the type and the kinematic regime of 
the processes which enter PDF fits. Therefore, their omission can lead to a significant 
distortion of the PDFs, thereby reducing their theoretical accuracy (see Ref. [25] for a 
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detailed discussion). In order to avoid this problem, it is necessary to supplement fixed- 
order calculations with all-order resummations of these large logarithms. 

Logarithmic enhancements of higher-order perturbative contributions originate from 
a number of different kinematic regions and require, in general, different resummation 
techniques (see e.g. Ref. [26] for a recent review). For instance, enhancements may take place 
when the centre-of-mass energy of the partonic collision is much higher than the hard scale 
of the process: this corresponds to the small-x region of the PDFs, and the resummation 
of such terms is known as high-energy or small-x resummation, see e.g. [27-30]. Small-x 
resummation is certainly relevant for PDF determination and might be needed to describe 
the most recent HERA data, where some tensions with fixed-order DGLAP have been 
reported [31-33]. A study of small-x resummation in PDF fits will be presented elsewhere. 
In this paper we concentrate instead on another type of logarithmic enhancement of higher 
order perturbative contributions which appear close to threshold for the production of the 
final states: this is the large-x kinematic region, and the resummation of logarithms from 
this region is known as large-x, soft gluon, or threshold resummation. 

All-order threshold resummations exist for many of the processes which play a central 
role in the exploration of the electroweak scale being pursued at the LHC. For instance, 
the current Higgs Cross Section Working Croup recommendation for the gluon-fusion cross 
section includes threshold resummation [34, 35], and resummed Higgs cross sections in this 
channel are available up to N^LL [36-38]. Additional resummed calculations for Higgs 
physics exist, for example for gluon-induced Higgs Strahlung [39] and for Higgs-pair pro¬ 
duction [40]. The basis for threshold resummation in QCD were laid in Refs. [41-43], in 
which explicit resummed expressions to NLL accuracy were given for processes relevant 
for global fits of parton distributions, such as deep-inelastic structure functions and Drell- 
Yan total cross section. Subsequently, higher-order resummed calculations have been made 
available for deep-inelastic scattering structure functions [44-48], invariant mass distribu¬ 
tions [38, 47-51] and rapidity distributions [52-58] in Drell-Yan production, and top quark 
pair production, both inclusive [59, 60] and differentially [61-63]. As far as processes rele¬ 
vant for New Physics searches are concerned, resummed calculations exist for squark and 
gluino production [64-67], stop quark pair production [68, 69], slepton and gaugino pair 
production [70-73] among many others. 

Moreover, it is well known [74, 75] that in the commonly used MS scheme, thresh¬ 
old resummation affects only partonic coefficient functions, while the singular part of the 
DCLAP splitting function is given, to any order in perturbation theory, by the cusp contri¬ 
bution: P{x,as) ~ rcusp(as)/(l — x), as X —)> 1. Therefore, to perform a resummed PDF 
fit it is only necessary to modify the partonic cross sections, while the NLO [76-79] and 
NNLO [80, 81] DCLAP evolution kernels remain unchanged. 

Despite the wide range of resummed calculations available, a complete global PDF fit 
including the effects of threshold resummation has never been produced (although some 
preliminary results were presented in Ref. [82]). A first study, restricted to non-singlet 
DIS structure functions, was performed in Ref. [83], finding that at NLO resummation 
could suppress the large-x valence quark PDFs by as much as ten percent. The impact 
of threshold resummation in direct photon production and its implications on the large-x 
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gluon was studied in Ref. [84], More recently, threshold resummation has been studied in 
the context of the CJ fits [85], with emphasis on the description of the large-x JLAB data. 
Given the impressive theoretical developments in the resummation of hard-scattering cross 
sections, it is clear that a state-of-the-art resummed global PDF fit is most timely. This is 
what we plan to achieve in this paper. 

To this end we will produce for the first time NLO+NLL and NNLO+NNLL threshold- 
resummed fits based on the NNPDF methodology [86-93]. Since NNPDF fitting is free from 
theoretical bias, due to the very flexible PDF parametrisation, it is sufficiently precise to 
be able to detect even small changes in PDFs due to threshold resummation of the various 
processes that go into the global fit. 

A major obstacle to producing a truly global resummed fit is that for a number of im¬ 
portant processes, in particular inclusive jet production and W production at the leptonic 
level, threshold resummation is not readily available. For inclusive jets, resummed calcu¬ 
lations have been used to determine approximate expressions [94-96] for the yet unknown 
NNLO contributions, but codes that provide all-order results are not publicly available. For 
W production, resummation is available only at the level of reconstructed W but not for 
the measured lepton-level distributions. 

For this reason, in this work we have begun by producing variants of the NNPDF3.0 
global NLO and NNLO fits [97] based only on those processes which can be consistently 
resummed: fixed-target and collider neutral and charged current deep-inelastic structure 
functions, fixed-target and collider neutral current Drell-Yan production, and inclusive top- 
quark pair production. These DIS+DY+top fits then provide a suitable baseline to compare 
with the NLL and NNLL resummed fits. One important drawback is that the resulting fits 
will be affected by larger PDF uncertainties as compared to the global NNPDF3.0 set, due 
to the missing experiments, affecting in particular gluon-initiated processes. In this respect, 
it will be important to produce updated resummed fits as soon as the missing resummed 
calculations become available. 

An important goal of this paper is to quantify the inaccuracies that affect current 
resummed calculations due to the inconsistent use of a fixed-order PDF with resummed 
partonic cross sections. As we will show, for final states with large invariant mass, close to 
the hadronic threshold, the main effect of the resummed PDFs is to bring the resummed 
hadronic calculation closer to the fixed-order result, thereby canceling partially the effect 
of the resummation in the matrix elements. On the other hand, for final states far below 
threshold, such as inclusive Higgs production at the LHC, the effect of the resummation 
on the PDFs can be small compared with the resummation in the matrix elements. We 
also find that, unsurprisingly, resummed and unresummed PDFs are much closer at NNLO 
than at NLO. Our results emphasise the need for a consistent use of resummed PDFs in 
resummed calculations: the use of fixed-order PDFs with resummed matrix elements can 
lead to misleading results, particularly at NLO. 

The outline of this paper is as follows. In Sect. 2 we review some basic concepts and 
results in threshold resummation, as well as their implementation. In Sect. 3, we discuss 
the settings of the global PDF fit used here to include threshold resummation effects, which 
is a variant of the recent NNPDF3.0 global fit. The results of the resummed fits are then 
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discussed in Sect. 4, where we compare resummed with fixed-order PDFs at NLO(+NLL) 
and NNLO(+NNLL). Then in Sect. 5 we discuss the implications of the resummed PDFs 
for LHC phenomenology, with emphasis on the mismatch that can arise if fixed-order PDFs 
are used in resummed calculations. Finally we summarise in Sect. 6 and discuss the delivery 
of the resummed PDFs produced in this work. 


2 Threshold resummation 


In this section we review the theoretical formalism of threshold resummation, and then we 
discuss its practical implementation in order to be able to use it in the resummed NNPDF 
fits. We work in the traditional framework of perturbative QCD (see e.g. [98]); alternative 
results can be obtained using the methods of Soft-Collinear Effective Theory (see Ref. [99] 
for a recent review and Refs. [100-103] for more detailed comparisons between the two 
formalisms). 


2.1 Theoretical framework 

We start by considering a hadron-level cross section 

a(x,Q‘^) = xY^ 

a^b 


f (- A-- ( 9l 

I ^ab ( ? /^F ) ^ab { ■> 2 ’ 2 I ’ 

Jx ^ / Z \ /ip / 


( 2 . 1 ) 


where a, b run over parton flavors, is the hard scale of the process, x is a dimensionless 
variable and x —>■ 1 defines the threshold limit. For the resummed fit we are going to 
consider three processes: deep-inelastic scattering (DIS) of a lepton off a hadron, the Drell- 
Yan process (DY) and top-anti-top production (tt). In DIS, is the off-shellness of the 
exchanged boson = —q^ and x = where p is the hadron momentum. In DY, Q is the 

zp-q ^ 

invariant mass of the lepton pair and x = ^, being ^/s the collider centre-of-mass energy. 
Finally, for ti, = 4mj and x = In Eq. (2.1), Cab{z,lA) is a parton luminosity, 
defined as 


^ab{z, /i ) 



fa fb{w,p^) , 

\W ) 


( 2 . 2 ) 


in the hadron-hadron collision case, while in the case of DIS it is just a single PDF. In the 
following we are going to set = fip = Q. 

In order to diagonalise the convolution integral, we take Mellin moments of Eq. (2.1): 


a{N,Q‘^)= / dxx 


N-2 
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where a 


In Mellin space the threshold limit corresponds to iV —)• oo and the aim of threshold 
resummation is to obtain a more reliable estimate of the hadron-level cross section by 
resumming to all orders in the strong coupling ag the logarithmically enhanced contributions 
to the partonic cross section aah at large N. The resummed partonic cross section can be 
written as the product of a Born contribution and an all-order coefficient function: 


cr 
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ab 
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where 


= '^90ab («^) exp5W(A^, a^), 

I 

{N, as) = In Aa -F In Ab + In Jc + In In (2.6) 


The notation ab —>■ cd has been chosen to accommodate all the processes that enter our 
fit. For tt production, we have to consider the resummation of two Born-level processes, 
namely qq —)• tt and gg —)• tt. For DIS instead we have V*q ^ q and for DY qq ^ V*. 
Moreover, while in DIS and DY we have one color structure, in the tt case we have two 
contributions, i.e. I = singlet, octet. 

Let us now examine the different contributions to the resummed exponent. If i is a 
color-singlet, then Aj = Jj = 1. For each initial-state QCD parton, we have an initial-state 
jet function 


f , ^ ^ ^ A / l 2\\ 

lnAi= / dz— - / -^Ai[as[q )), i = a,b. 

Jo ^ - z q 

For each massless final-state QCD parton we have a final-state jet function 
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( 2 . 8 ) 

while there is no jet-function for t or t. Finally we also have a large-angle soft contribution, 
which depends in principle on both the process and the color flow: 


^^^aLcd= I dz^ 




N-l 


1 - Z 


— D™ 


ab^cd 


[as mi-zf)). 


(2.9) 


The functions Ai (a^), Bi (a^), D® (a^), and obviously ^o®(Q(s); are free of large loga¬ 
rithms and can be computed in fixed-order perturbation theory. The accuracy of their 
determination fixes the logarithmic accuracy of the resummation. In particular, (N)NLL 
requires Ai to second (third) order in the strong coupling a^, and Bi, Df'\ and to 
first (second) order.^ Threshold resummation is actually known to N^LL for DIS [45] and 
DY [38, 49-51], and to NNLL for ti production [59, 104]. 

^This accuracy is sometimes referred to as (N)NLL'. For a precise definition of all possible accuracies 
and their nomenclature in threshold resummation, see Tab. 1 of Ref. [36]. 
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We have left out of our discussion inclusive jet production.^ The general framework to 
perform this resummation has been worked out long ago [105]. However, different treat¬ 
ments of the jet kinematics at threshold can lead to substantially different results, see e.g. 
Ref. [94] and [106]. Moreover, depending on the way the threshold limit is defined, NLL 
resummation can be affected by non-global logarithms [107] and the result may acquire a 
non-trivial dependence on the jet algorithm [108, 109]. In addition, as previously mentioned, 
computer programs that implement threshold resummation for jet production are not, to 
the best of our knowledge, publicly available. On the other hand, recent progress [110-112] 
has shown that NNLL accuracy is perhaps achievable in the near future. We leave a detailed 
phenomenological analysis of jet production and its inclusion in a PDF fit to future work. 

The Mellin integrals in the resummed expression Eq. (2.6) are often evaluated in the 
N ^ oo limit, thereby keeping only those contributions that do not vanish at large N 
and behave as powers of InN. In this approximation, which we refer to as Wsoft in the 
following, the resummed coefficient function becomes 


a^) = go{as)ex.pS{\nN,as), 


5(lnA^, 


—fi(i(aslnA^) + g 2 {as InA^) -h 0 ^ 5 ( 3(05 In A^) -k ... 

(y,g 


( 2 . 10 ) 


where, in order to simplify our notation, we henceforth drop all the flavor and color-flow 
indices and it is understood that all the modifications we discuss are applied to each partonic 
subprocess and each color-flow. The functions gi^aglnN) with i > 1 resum a”ln"'A^ 
contributions to all orders in perturbation theory. They can be derived directly by the 
integral representations Eqs. (2.7)-(2.9) by computing the integrals as an expansion in 
powers of as at fixed in the large N limit. 

Other resummation schemes, which are equivalent to Wsoft at large N, but preserve the 
analytic structure of fixed-order coefficient functions at finite N, have been considered in the 
context of Higgs production [36, 113, 114] and heavy quark production [115]. The extension 
of these resummation schemes to DIS and DY, and their application to the determination 
of PDFs, will be considered elsewhere. 

Threshold resummation can be extended to rapidity distributions, see e.g. [52-58]. In 
this work we follow the approach of Ref. [58]. The basic observation is that the resummed 
partonic rapidity distribution coincides with the rapidity integrated one up to terms which 
are power-suppressed in the threshold limit. Therefore, in order to obtain the hadron-level 
resummed rapidity distribution, we have only to modify the parton luminosity. While we 
refer the Reader to Ref. [58] for an explicit derivation, here we limit ourselves to note that 
the resummed rapidity distribution is constructed in such a way that the integral over 
rapidity gives back the resummation of the rapidity-integrated cross section. 

Finally, we mention that the calculation of hadron-level cross sections and distribu¬ 
tions from resummed results in N space requires a prescription because of the presence of 
a logarithmic branch-cut for real N > Nl, originating from the Landau pole of the running 
coupling. As a consequence, the resummed result does not admit an inverse-Mellin trans¬ 
form. Different solutions to this problem exists, such as the Minimal Prescription [116], 

^We acknowledge discussions with Mrinal Dasgupta and Werner Vogelsang on this topic. 
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which consists on a simple modihcation of the Mellin inversion integral, and the Borel pre¬ 
scription [58, 117-120], which relies on a Borel summation of the divergent series of the 
order-by-order inverse Mellin transform of the resummed coefficient function. In this paper 
we adopt the Minimal Prescription, but we stress that from a practical point of view, dif¬ 
ferences between these prescriptions become only relevant at extremely large values x [120], 
a region where no experimental data is available. 

2.2 Numerical implementation 

In this section we discuss the numerical implementation of the N-soft threshold resum¬ 
mation described above. For the PDF hts performed in this work, the processes that we 
are interested in are DIS (both neutral and charged currents), lepton-pair invariant mass 
and rapidity distribution for Drell-Yan production, and inclusive top pair production cross 
section. For DIS and Drell-Yan, we use a new version of the public code ResHiggs, written 
originally [36] to perform threshold resummation of Higgs inclusive cross section, including 
several improvements with respect to standard N-soft resummation, and later extended [37] 
to also perform (improved) resummation in the Soft-Collinear Effective Theory formalism. 
Because of the inclusion of additional processes, the new version of this code changes name 
from ResHiggs to TROLL, standing for TROLL Resums Only Large-x Logarithms, publicly 
available at the webpage [121]. To give continuity with the original code ResHiggs, the hrst 
version of TROLL is v3.0. For top pair production we use the public code Top++ [7]. 

The code TROLL is designed to provide only the contribution of the resummation, 
while the hxed-order calculation is obtained from a separate code (in our case, the same 
FKgenerator code used for the NNPDF3.0 hts). More specihcally, the output of TROLL is 
dehned as the difference between a resummed iF-factor at N-^LO+N^LL and a 
hxed-order 7C-factor at N-^LO, such that 

'^NYO+N'^LL = '^NYO + X ^J-^N^LL > (2-11) 

where all the cross sections appearing in the above equation are evaluated with a common 
N-^LO+N^LL PDF set. Internally, is computed by subtracting off the expansion 

of the resummed coefficient up to 0(ai) from the coefficient itself, multiplying this by 
the parton luminosity, computing the inverse Mellin transform and hnally dividing by the 
parton luminosity. 

There are several advantages in using the Aii'-factors dehned in Eq. (2.11): the hxed- 
order normalisation is irrelevant, the iF-factor is much less sensitive to the input PDFs than 
the cross section itself, and hnally, since the resummed contribution has the same kinematic 
structure as the Born cross section (soft radiation does not change the kinematics), the effect 
of phase space constraints like kinematic cuts are correctly taken into account if they are 
applied to the LO cross section in Eq. (2.11). 

We note that in DIS, Target Mass Corrections (TMCs) at next-to-leading twist are 
included in the resummation according to the same prescription used in the NNPDF htting 
code [88], which amounts to multiplying the Mellin transforms of the partonic coefficient 
functions by an N-dependent factor. No TMCs are included for the hxed-target DY data. 
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Neutral current DIS 


Neutral current DIS 




Figure 1. Aif-factors Eq. (2.11) for the neutral current DIS structure function F 2 {x,Q), as a 
function of a:, for Q = 2 GeV and Q = 30 GeV. The plot on the left corresponds to j = 1, fc = 1 in 
Eq. (2.11), i.e. NLO and NLL, while the one on the right to j = 2, k = 2, i.e. NNLO and NNLL. 
The effect of adding TMGs is shown as a thin solid line. 



Figure 2. Same as Fig. 1 for the NLL (left plot) and NNLL (right plot) resummation of neutral- 
current Drell-Yan invariant mass distribution at the Tevatron and at the LHG. 


In Fig. 1 we show the AiF-factors for the neutral current DIS structure function 
F 2 {x,Q), as a function of x, for Q = 2 GeV and Q = 30 GeV. The plot on the left 
corresponds to j = 1, A: = 1 in Eq. (2.11), i.e. NLL to be matched to NLO, while the 
one on the right to j = 2, A: = 2, i.e. NNLL to be matched to NNLO. We note that the 
resummation enhances the cross section at large x, while it gives a very small contribution 
at small x, as it should. We also note a dip in the region of intermediate x, which is also 
present in hxed-order calculations [122]. 

TMG effects are also shown as light shadows to the actual curves: as expected, at large 
scales they are negligible, while at smaller scales their effect is non-negligible, in particular 
at large x, where they reduce the effect of the resummation. Note that in the dehnition of 
AK we use the fact that the same TMGs are already included in the LO cross section, so 
much of their effect cancels out. 

In Fig. 2 we show the corresponding AiL-factors this time for Drell-Yan invariant mass 
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Figure 3. Same as Fig. 1 for the neutral-current Drell-Yan rapidity distribution, for different 
experiments and different values of the lepton invariant mass. 



LHC 7 TeV 

LHC 8 TeV 

O'NLO-I-NLl/o'NLO 

1.086 

1.081 

CTNNLO-l-NNLL/crNNLO 

1.031 

1.029 

crNNLo/<7NLO 

1.123 

1.122 


Table 1. iF-factors for tt production at LHC at 7 and 8 TeV. The first line is obtained using NLO 
PDFs, while the second and third lines are obtained with NNLO PDFs (in all cases, both numerator 
and denominator are computed with the same PDFs). 


distributions, for LHC at -y/s = 7 TeV and for the Tevatron at -^/s = 1.96 TeV. In Fig. 3 
we also present the results for the lepton pair rapidity distribution in DY processes as a 
function of Y with Vnax = \ Y^{s/M‘^) being the maximum rapidity of the lepton 
pairs allowed by kinematics. The experiments cover different kinematic regimes: close to 
threshold (Fermilab’s fixed-target Drell-Yan experiments at M = 7 GeV), an intermediate 
region (Tevatron and LHC at the Z pole) and away from threshold (low mass DY at LHC). 
We note that threshold resummation always gives a significant correction at large rapidities. 
For the hxed-target kinematics and M = 7 GeV, the effect of threshold resummation is 
important even at central rapidities. 

Finally, in Table 1 we collect the iC-factors for tt production obtained using Top++. In 
this case we provide directly the iC-factors for the (N)NLO+(N)NLL over (N)NLO cross 
sections, with both numerator and denominator computed with the same (N)NLO PDFs. 
We can see from the table that the impact of the resummation is non-negligible, and in 
fact rather important especially at NLO+NLL, where the correction is about 9% of the 
hxed-order NLO result, comparable to the NNLO correction. Even at NNLO+NNLL, the 
effect of resummation is comparable to other theory uncertainties like the values of as(m^) 
or of the top quark mass [123]. 
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3 Settings of the resummed PDF fit 


In this section we present the settings used in the resummed PDF hts. These are con¬ 
structed as variant of the recent NNPDF3.0 global hts [97]: they use exactly the same 
htting methodology, the same input parameters (strong coupling, heavy quark masses, 
etc.), and the same hxed-order theoretical calculations. The experimental dataset is also 
similar except that some specihc processes have been excluded. 

In this section, we hrst review the experimental data that can be consistently included 
in a threshold resummed global PDF analysis, and then we explain the procedure used 
to construct the resummed iP-factors that will be used to include threshold resummation 
in the NNPDF hts. We also show the resulting resummed iP-factors for a representative 
subset of the experiments used in the ht. 

3.1 Experimental data 

In a PDF ht with threshold resummation, as compared to hxed-order hts, some datasets 
cannot be included since for these processes threshold resummation is either unknown or 
not currently available in a format that can be used in a ht. In particular, when compared 
to NNPDF3.0, in the present resummed ht we include all the neutral and charged current 
DIS data, neutral current DY production and top quark pair production data. However we 
exclude the DY charged current datasets, for which data is provided in terms of the lepton 
kinematics, and the construction of resummed expressions is more involved, and inclusive 
jet production for the reasons discussed in Sect. 2.1. 

In Table 2 we list all the datasets used in the NNPDF3.0 NLO and NNLO global anal¬ 
ysis, and indicate whether or not they are now included in the NLL and NNLL resummed 
NNPDF3.0 hts. For each dataset we also display the corresponding measured observable, 
and the relevant publication. A more complete description of each of these datasets, as well 
as of their impact in terms of PDF constraints in the global ht, can be found in [97]. 

From Table 2, we infer that, when compared to the global ht, the resummed hts lose 
experimental constraints on the medium and large-x gluon (due to the exclusion of the jet 
data) and on the quark-havor separation (due to the exclusion of the W data). Still, given 
that we include in the resummed ht more than 3000 data points, the loss of accuracy due to 
the exclusion of these datasets is not dramatic, as we will show in Sect. 4. In future studies, 
we aim to include the missing processes once the corresponding resummed calculations 
become available. 

The kinematic cuts applied in the present hts closely follow the ones of the NNPDF3.0 
hxed-order analysis. In particular, a cut on the hnal-state invariant mass of DIS data 

> 12.5 GeV^ is applied, in order to reduce the dependence on higher-twists at large x. 
It would be interesting to loosen this cut in future analyses, in order to test the stability of 
the leading-twist PDF determination once the large-x resummation is included; this might 
also allow us to include additional large-x, low-Q^ DIS measurements, such as for example 
JLAB data [164]. 

In addition, a stability analysis of the calculations for neutral-current Drell-Yan pro¬ 
duction indicates that our results, even when supplemented with resummation, become 
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Experiment 

Observable 

Ref. 

NNPDF3.0 global 

NNPDF3.0 DIS+DY+top 




(N)NLO 

(N)NLO ]+(N)NLL] 

NMC 


[124, 125] 

Yes 

Yes 

BCDMS 

Fi,Fi 

[126, 127] 

Yes 

Yes 

SLAG 

Fi,F^ 

[128] 

Yes 

Yes 

CHORUS 


[129] 

Yes 

Yes 

NuTeV 

CC, charm 
'^vN 

[130] 

Yes 

Yes 

HERA-I 

NC CC 
•^dis 1 ^dis 

[131] 

Yes 

Yes 

ZEUS HERA-II 

NC CC 

^dis 5 ^dis 

[132-135] 

Yes 

Yes 

HI HERA-H 

^NC ^CC 
^dis 1 ^dis 

[136, 137] 

Yes 

Yes 

HERA charm 

NC, charm 
'^dis 

[138] 

Yes 

Yes 

DY E866 

^NC /^NC 

^DY,p? ^DY,d/*^DY,p 

[139-141] 

Yes 

Yes 

DY E605 

_NC 

'^DY,p 

[142] 

Yes 

Yes 

CDF Z rap 

_NC 

'^DY,p 

[143] 

Yes 

Yes 

CDF Run-H kt jets 

O'jet 

[144] 

Yes 

No 

DO Z rap 

NC 

^DY,p 

[145] 

Yes 

Yes 

ATLAS ^ 2010 

'^DY,p 

[146] 

Yes 

Yes 

ATLAS W 2010 

'^DY,p 

[146] 

Yes 

No 

ATLAS 7 TeV jets 2010 

Ojet 

[147] 

Yes 

No 

ATLAS 2.76 TeV jets 

Ojet 

[148] 

Yes 

No 

ATLAS high-mass DY 

NC 

^DY,p 

[149] 

Yes 

Yes 

ATLAS W pT 

^DY,p 

[150] 

Yes 

No 

CMS W electron asy 

^cc 

'^DY,p 

[151] 

Yes 

No 

CMS W muon asy 

'^DY,p 

[152] 

Yes 

No 

CMS jets 2011 

Ojet 

[153] 

Yes 

No 

CMS W + c total 

NC, charm 
'^DY.p 

[154] 

Yes 

No 

CMS 2D DY 2011 

^NC 

'^DY,p 

[155] 

Yes 

Yes 

LHCb W rapidity 

'^DY,p 

[156] 

Yes 

No 

LHCb Z rapidity 

^NC 

'^DY,p 

[157] 

Yes 

Yes 

ATLAS CMS top prod 

a{tt) 

[158-163] 

Yes 

Yes 


Table 2. List of all the experiments that were used in the NNPDF3.0 global analysis, and whether 
or not they are now included in the present (N)NLL resummed fits (and in the corresponding 
baseline fixed-order fits). For each dataset we also provide the type of cross section that has been 
measured and the corresponding publication references. 


unstable for data points too close to the production threshold, either because the invari¬ 
ant mass Mil of the Drell-Yan pairs is too large, or because the rapidity Y is too close to 
the kinematic boundary Ymax- Therefore, we have supplemented the NNPDF3.0 kinematic 
cuts with two additional cuts for the fixed-target Drell-Yan experiments, as summarised in 
Table 3. For the collider Drell-Yan data, the cuts are the same as in NNPDF3.0. 

It is useful to quantify which experiments determine the behaviour of the large-x PDFs 
in the global analysis. In Fig. 4 we compare the relative PDF uncertainties in variants of 
the NNPDF3.0 NLO fit based on different input datasets: HERA-only, no LHC data, no 
jet data, a DIS+DY+top fit (our baseline for the resummed fits) and the global fit. We 
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Experiment 

Kinematic cuts 

DIS 

> QLn = 3.5 GeV2 
VF2 > = 12.5 GeV2 

Fixed target Drell-Yan 

T < 0.08 

|Y|/Y„,ax < 0.663 


Table 3. Kinematic cuts applied to the DIS and fixed-target Drell-Yan data in the baseline and 
resummed fits. For the collider Drell-Yan data, the cuts are the same as in NNPDF3.0. 


NNPDF3 NLO, 0^=10'* GeV^ 


NNPDF3 NLO, 0^=10“ GeV=^ 



“b 


0.8 

0 


- HERA-only 

- No LHC data 

.No jet data 

. DIS+DY+top 

-Global 


NNPDF3 NLO, 0^=10“ GeV^ 


NNPDF3 NLO, 0^=1 O'* GeV= 




Figure 4. Comparison of the relative PDF uncertainty at large-a; between NNPDF3.0 NLO fits 
based on different input datasets. We show the results for the gluon, total quark singlet, total 
valence and u quark PDFs, at a typical LHC scale of = 10"^ GeV^. The fits shown are the 
HERA-only, no LHC data, no jet data, DIS+DY+top (our baseline for the resummed fits) and 
finally the global fit. 


focus in the large-x region for = 10^ GeV, a typical scale for LHC phenomenology. We 
show the gluon, the total quark singlet, the total valence and the u quark PDFs. From 
this comparison we see that the PDF that at large-x is most dependent on the choice of 
input dataset is the gluon. For the total valence and singlet quark PDFs, the bulk of the 
constraints are provided by the DIS and fixed target Drell-Yan data, which are common 
datasets in all these fits (except for the HERA-only fit). For the sea quarks, in this case the 
u quark, the information on both jet data and LHC data are necessary to achieve the best 
possible accuracy. The baseline DIS+DY+top fit is slightly less accurate at large-x for the 
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quark flavor separation due to the missing charged-current Drell-Yan data. 


3.2 Calculation of resummed K-factors 

As mentioned in the introduction, in the MS scheme all the effects of threshold resummation 
are encoded in the partonic cross sections, and thus parton evolution is the same as in 
fixed-order calculations. Therefore, apart from the modification of the hard-scattering 
cross sections, all theoretical settings in the resummed fit will be the same as those of 
the NNPDF3.0 fixed-order analysis, including the use of the FONLL general-mass VFN 
scheme [165], the values of the heavy quark masses, and so on. We will produce results for 
single a value of the strong coupling, as{rn?z) ~ 0.118. 

As discussed in Ref. [97], the NNPDF3.0 global analysis for hadronic observables always 
uses fast NLO calculations [166-168] supplemented with NNLO/NLO iF-factors when re¬ 
quired. These are defined as the ratio of the NNLO over the NLO bin-by-bin cross sections, 
using a common PDF luminosity computed with a NNLO PDF set. For the resummed fits, 
we follow exactly the same procedure: we include the effect of resummation supplementing 
the fixed-order computation with a AT-factor. Since the iL-factor is computed externally 
using a fixed set of PDFs, the fit is re-iterated several times, recomputing each time the 
resummed AT-factor using as input PDFs those obtained from the previous iteration. 

As discussed in Sect. 2.2, the resummed contributions for DIS and DY processes is 
obtained using the program TROLL in the form of AAT-factors, Eq. (2.11), and hence must 
be converted into actual AT-factors. For DIS cross sections, since the NNLO calculation is 
implemented exactly in the NNPDF fitting code, this is done according to 


^N*LO+N'=LL 
7^N'“L0+N'=LL _ ^ 

^DIS - ^n^LO 


1 + 


a 


LO 


a 


N'^LO ’ 


(3.1) 


with A: = 1,2 for NLO+NLL and NNLO+NNLL respectively. For hadronic processes we 
use a similar expression, but (at NNLO) also including the NNLO/NLO A'-factor, 


^NLO+NLL _ 
-^hadr — 


7^NNL0-I-NNLL _ 
-^hadr — 


a 


NLO+NLL 


NLO 


a 


a 

NNLO+NNLL 


= 1-1- AiATnll • 


a 


LO 


NLO 


a 


a 


NLO 


= K 


NNLO 


+ A2A'nNLL 


a 


LO 


NLO 


a 


(3.2) 

(3.3) 


where = ^t^nlo^'^nlo above expressions, all contributions are meant to 

be computed with the same N^LO+N^LL PDF set. For the leading-order cross section 
cr'"®, dedicated FK tables with LO coefficient functions but NLO and NNLO PDF evolution 
have been produced using the APFEL program [169], validated with the same FKgenerator 
internal code used in the NNPDF3.0 fits. 

In principle, all these contributions should be recomputed at each iteration of the fit; 
in practice, the computation of /f^NLO jg consuming, so for this contribution we use a 
fixed value. Specifically, in the present work, these fixed-order NNLO/NLO LC-factors are 
the same as in the NNPDF3.0 fits, with the exception of those for fixed-target Drell-Yan 
experiments, which have been recomputed using Vrap [170] with the NNPDF3.0 global PDF 
set as input. 
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Figure 5. The resummed X-factors for DIS, Eq. (3.1), for a representative subset of the experiments 
included in the resummed fit: BCDMS NMC cr^c. ^LAC and CHORUS a^N- We show 
both the results corresponding to NLL and to NNLL resummation. The DIS kinematics {x,Q^,y) 
is that of the corresponding experimental data included in the fit, plotted just as a function of x, 
so for each value of x there are measurements at different values of and y. 


For the computation of the resummed iC-factors Eq. (3.1)-(3.3), we find that two 
iterations of the fit are enough to reach a satisfactory convergence, meaning that these 
iC-factors are essentially unchanged if we use resummed PDFs from the last or from the 
next-to-last iteration of the fit. 

It is now interesting to illustrate the effect of the (N)NLL resummation for some of the 
datasets used in the present resummed PDF fit. To this purpose, we plot the DIS, Eq. (3.1), 
and hadronic, Eqs. (3.2) and (3.3), resummed iC-factors for representative experimental 
datasets with exactly the same kinematics as for the data points to be used in the fit. In 
these calculations, we have consistently used the NNPDF3.0 DIS+DY+top NLO+NLL and 
NNLO+NNLL PDF sets discussed in the next Section, with as{rn?z) = 0.118, in both the 
fixed-order and resummed cross sections. To isolate the effect of the resummation, in these 
comparison plots we will factor out Rj^NLO hadronic NNLO+NNLL resummed 

iP-factor Eq. (3.3). Note that in these plots we will only include those data points that 
satisfy the kinematic cuts imposed in the fit, summarised in Table 3. 

The results for the DIS case are shown in Fig. 5. For each experiment, we show 
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both the results corresponding to NLL and to NNLL resummation. The DIS kinematics 
(x, y) is that of the associated experimental measurements, so for each value of x there 
are measurements at different values of and y. We do not show the results for any of 
the HERA datasets, for which the effect of the resummation turns out to be negligible since 
the data is either at small-x or at high scales. As expected, the impact of the resummation 
is only relevant at large x, and of course the impact of the resummation decreases when 
more fixed perturbative orders are included in the calculation. 

From the results of Fig. 5 we see that the effect of threshold resummation is most 
important for the BCDMS data, while it is milder for the other experiments. Effects of 
NLL resummation reach up to 15% at the highest values of x available, which is reduced 
to up to a few percent for NNLL resummation (since part of the effects at NLL are now 
included in the fixed-order NNLO calculation). Note also that the cut in removes most 
of the large-x SLAG data, where resummation effects are very large. We also note that the 
effect of resummation is comparable to the experimental uncertainties, and thus we should 
expect to see an impact on the resulting large-x parton distributions. 

The corresponding results for the resummed AT-factors Eq. (3.3) for representative 
neutral-current DY experiments are collected in Fig. 6. As in the case of DIS, except for Z 
peak measurements, for each DY rapidity Y value there are various data points at different 
invariant mass M. We show the Drell-Yan E866 pp cross sections, the CDF Z rapidity 
distribution, the CMS double-differential Drell-Yan distribution at 7 TeV and the LHCb 
Z —7- rapidity distribution. 

From Fig. 6 we verify the expectation that the impact of resummation is rather more 
important at NLO+NLL than at NNLO+NNLL, and that it grows with the di-lepton 
rapidity (since in this case the kinematic threshold is approached). For those collider 
measurements differential in rapidity, the effect of (N)NLL resummation can be as large 
as 50% (20%) at the highest rapidities. For the fixed-target DY experiments the effect 
of resummation is substantial even at NNLL, since in this case many data points have 
kinematics close to threshold. For example, for the E866 pp dataset, the effect of the 
resummation results is an enhancement of the cross section that can be as large as 35% at 
NLL, and 20% at NNLL. 

Following this discussion on the settings used to produce the threshold resummed fits, 
in the the next section we turn to explore the actual effects that the inclusion of resummed 
calculations have on the NNPDF3.0 PDFs. 


4 Parton distributions with threshold resummation 

In this section we discuss the results of the NNPDF3.0 fits with threshold resummation. 
One important difference of the resummed fits as compared to the NNPDF3.0 global fits is 
that the dataset is different, because we leave out the inclusive jet and W production data, 
as discussed in Sect. 3. Therefore, first we quantify the information loss due to the reduced 
dataset by comparing the global NNPDF3.0 fits and the reduced dataset fits obtained with 
fixed-order matrix elements (henceforth denoted as the baseline fits). 
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Figure 6. Same as Fig. 5 for selected Drell-Yan experimental datasets included in the fit: the E866 
pp neutral Drell-Yan cross sections, the CDF Z rapidity distribution, the CMS double-differential 
Drell-Yan distribution at 7 TeV and the LHCb Z pp rapidity distribution. The resummed K- 
factors are now those defined in Eqs. (3.2) and (3.3), but to isolate the effect of resummation from 
that of the fixed-order NNLO corrections, in the NNLL case we divide Eq. (3.3) by 
two left plots data points differ by the values of the rapidity and the invariant mass of the pair, but 
only the dependence on the rapidity is shown. 


Having established this, we move to quantify the impact of threshold resummation on 
the fit quality and the resulting PDFs, by comparing fits at NLO and NLO+NLL first, and 
then at NNLO and NNLO+NNLL. This is done both for DIS-only fits and for DIS+DY+top 
fits, and the comparison is performed both at the level of PDFs and of Finally, we assess 
the impact of threshold resummation at the level of partonic luminosities. 

The phenomenological implications of the resummed PDFs for LHC applications will 
be discussed in the next section. 

4.1 Baseline fixed-order fits 

First of ah, we compare the baseline fixed-order fits with the NNPDF3.0 global sets. As 
mentioned in Sect. 3, in the fixed-order baseline fits ah settings are identical to those of 
NNPDF3.0 with the only difference of the use of a reduced dataset, see Table 2. Therefore, 
we expect the two fits to be consistent, with the baseline fit affected by larger PDF uncer- 
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NNPDF3.0 NLO, 0^=10“ GeV^ 


NNPDF3.0 NLO, 0^=10“ GeV^ 



Figure 7. Comparison of the fixed-order NNPDF3.0 NLO fits based on different datasets: global, 
DIS-only and DIS+DY+top, for as{rn^) = 0.118, at a typical LHC scale of = 10“^ GeV^. Results 
are normalised to the central prediction of the NNPDF3.0 NLO global fit. From left to right and 
from top to bottom, we show the gluon, the total quark singlet, the total valence PDF and the total 
strangeness. 


tainties due to the reduced dataset. For simplicity, we restrict these comparisons to NLO, 
since the impact of the reduced dataset is roughly independent of the perturbative order. 

We have produced two baseline fits: one with the all the data marked in the last column 
of Table 2, and the other with only DIS-data included. In Fig. 7 we compare the NNPDF3.0 
NLO DIS-only and DIS+DY+top set with as{m‘^) = 0.118, with the corresponding global 
set.^ In both cases we use Wep = 100 replicas, and the comparison is performed at a typical 
LHC scale Q = 100 GeV. Results are normalised to the central prediction of the global fit. 
From left to right and from top to bottom, we show the gluon, the total quark singlet, the 
total valence PDF and the total strangeness. 

As we can see, there is a reasonable agreement for most of the PDF flavours and of 
momentum fraction x between the three fits, with as expected larger PDF uncertainties 
in the DIS-only and DIS+DY+top fits due to the reduced dataset. The DIS-only fit is 
relatively close to the global fit for the large-x quarks, since these are well constrained by 
the DIS fixed-target data. On the other hand, the DIS-only fit is affected by rather larger 

®In the rest of this section, we concentrate only on the large-® region of the PDFs, since as we will show 
the effects of threshold resummation are negligible at medium and small-®. 
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uncertainties as compared to the global fit for the gluon (due to the missing jet data) and 
for the total valence (due to the missing Drell-Yan data that constrains flavor separation). 
In any case, the DIS-only and the global ht are always consistent at the one-sigma level. 

Concerning the DIS+DY+top ht, for quark PDFs (singlet, valence and strangeness) 
the results of the DIS+DY+top ht are quite close to the global ht. Therefore, we can 
conclude that for quark-initiated processes, calculations done with the DIS+DY+top hts 
are essentially equivalent to those performed using the global PDFs. The only differences 
are as expected related to the gluon PDF, where the missing inclusive jet data cause a 
substantial increase in the large-x gluon PDF uncertainties compared to the global ht. 
Note however that in the resummed DIS+DY+top ht a handle on the large-x gluon is still 
provided by the total top-quark pair production cross section [123, 171]. 

After having established the impact of the reduced datasets on the baseline hts that 
will be used for the resummation, in the following we concentrate in quantifying the impact 
of resummation for hts based on a common dataset, hrst for the DIS-only hts, and then for 
the DIS+DY+top hts. 

4.2 DIS-only resummed PDFs 

Now we present the results of the resummed hts. We begin with the DIS-only hts, and 
compare the baseline NLO and NNLO hxed-order with the corresponding NLO+NLL and 
NNLO+NNLL threshold resummed hts. 


Experiment 

NNPDF3.0 DIS-only 


NLO 

NNLO 

NLO+NLL 

NNLO+NNLL 

NMC 

1.36 

1.39 

1.36 

1.32 

SLAG 

1.12 

1.15 

1.02 

1.04 

BCDMS 

1.19 

1.20 

1.21 

1.22 

CHORUS 

1.10 

1.05 

1.09 

1.07 

NuTeV 

0.52 

0.46 

0.55 

0.51 

HERA-I 

1.07 

1.13 

1.06 

1.07 

ZEUS HERA-II 

1.40 

1.42 

1.42 

1.43 

HI HERA-H 

1.67 

1.79 

1.68 

1.74 

HERA charm 

1.28 

1.29 

1.29 

1.24 

Total 

1.237 

1.257 

1.237 

1.242 


Table 4. The per data point for all experiments included in the DIS-only threshold resummed 
hts, at NLO and NNLO, compared with their resummed counterparts. 


In Table 4 we provide the per data point for all experiments included in the DIS- 
only threshold resummed hts, at NLO+NLL and NNLO+NNLL, to be compared with 
their unresummed counterparts.^ From Table 4 we see that, as expected, the impact of 
resummation is moderate and restricted to the hxed-target DIS experiments. In the case of 
SLAG, there is a clear improvement in the due to the inclusion of threshold resummation, 

^ As in NNPDF3.0, the present fits used the to definition of the covariance matrix for the minimisation, 
but then use the experimental definition to assess the consistency between theory and data. See [2, 172] 
for the explanation of the different definitions of the x^ estimators. 
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NNPDF3.0 DIS-only, Q^=10^ GeV^ 



NNPDF3.0 DIS-only, 0^=10'* GeV^ 



NNPDF3.0 DIS-only, 0^=10“ GeV^ 



NNPDF3.0 DIS-only, Q^=10^ GeV^ 



Figure 8. Comparison between the NNPDF3.0 DIS-only NLO fit and the corresponding NLO+NLL 
fit, for as{rn\) = 0.118, at a typical LHC scale of = 10^ GeV^. PDFs are normalised to the 
central value of the fixed-order fit. 


both at NLO and at NNLO. For other experiments the change in is not significant, 
meaning that the small effect of threshold resummation can be absorbed in the fitted 
PDFs. The total is slightly improved when going from the NNLO to the NNLO+NNLL 
fit, while it is essentially unaffected in the NLO+NLL case. 

Now we turn to study the impact of the resummation on the PDFs themselves. In 
Fig. 8 (for the NLO) and in Fig. 9 (for the NNLO) we compare the NNPDF3.0 DIS-only 
(N)NLO set with Qis(m^) = 0.118, with the corresponding (N)NLO+(N)NLL threshold 
resummed PDFs, respectively. 

First of all, we note that as expected the inclusion of threshold resummation affects only 
PDFs at large x, for x > 0.1, which is consistent with the modifications that resummation 
induces on the DIS structure functions. We also see that the impact on the PDFs is more 
important at NLL than at NNLL, again as expected since NLL captures part of the NNLO 
corrections to the DIS structure functions. For the quark PDFs, the effect of resummation 
is a suppression of the central values for quite large x. For example, for the NLO+NLL 
fit, the total quark singlet S(x,(5^) is suppressed by ~ 5% at x ~ 0.6. One also observes 
a small enhancement of the valence PDF for x ~ 0.2, presumably due to a compensation 
for the suppression at very large x through the valence sum rules. Therefore, we expect 
resummation to have phenomenological impact for the calculation of quark-initiated heavy 
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NNPDF3.0 DIS-only, Q^=10^ GeV^ 



NNPDF3.0 DIS-only, 0^=10'* GeV^ 



NNPDF3.0 DIS-only, 0^=10“ GeV^ 



NNPDF3.0 DIS-only, Q^=10^ GeV^ 



Figure 9. Same as Fig. 8 now comparing the NNLO DIS-only fit with the corresponding 
NNLO+NNLL fit. 


production processes in BSM scenarios, which probe rather large values of x. The large-x 
gluon is also suppressed at large-x, though PDF uncertainties are too large in a DIS-only 
fit to make this suppression significant. 

4.3 DIS+DY+top resummed PDFs 

Now we turn our discussion to the impact of threshold resummation in the case of the fits 
based on the DIS+DY+top dataset. In Table 5 we provide the again using the exper¬ 
imental definition, for all experiments included in the DIS+DY+top threshold resummed 
fits, at NLO+NLL and NNLO+NNLL, compared with their fixed-order counterparts. 

As we can see from Table 5, in the NLO+NLL fits the fit quality of most of the 
experiments is improved as compared to the fixed-order NLO fit. This is especially marked 
in the case of SLAG, as in the DIS-only fit, but also for the CHORUS neutrino structure 
functions, the CDF Z rapidity distribution, ATLAS high-mass DY, LHCb Z rapidity and 
the top quark pair production. The only exception is the fixed target Drell-Yan data, 
where resummation makes the worse. Note however that in the resummed fit the overall 
balance between experiments in the global fit is modified compared to the fixed-order fit, 
so this does not necessarily imply that resummation degrades the internal fit quality for 
this particular observable.^ At the level of total see that fixed-order and resummed 

®We have checked that in a fit based only on HERA data and fixed-target Drell-Yan data, in both the 
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Experiment 

NNPDF3.0 DIS+DY+top 


NLO 

NNLO 

NLO+NLL 

NNLO+NNLL 

NMC 

1.39 

1.34 

1.36 

1.30 

SLAG 

1.17 

0.91 

1.02 

0.92 

BCDMS 

1.20 

1.25 

1.23 

1.28 

CHORUS 

1.13 

1.11 

1.10 

1.09 

NuTeV 

0.52 

0.52 

0.54 

0.44 

HERA-I 

1.05 

1.06 

1.06 

1.06 

ZEUS HERA-II 

1.42 

1.46 

1.45 

1.48 

HI HERA-H 

1.70 

1.79 

1.70 

1.78 

HERA charm 

1.26 

1.28 

1.30 

1.28 

DY E866 

1.08 

1.39 

1.68 

1.68 

DY E605 

0.92 

1.14 

1.12 

1.21 

CDF Z rap 

1.21 

1.38 

1.10 

1.33 

DO Z rap 

0.57 

0.62 

0.67 

0.66 

ATLAS ^2010 

0.98 

1.21 

1.02 

1.28 

ATLAS high-mass DY 

1.85 

1.27 

1.59 

1.21 

CMS 2D DY 2011 

1.22 

1.39 

1.22 

1.41 

LHCb Z rapidity 

0.83 

1.30 

0.51 

1.25 

ATLAS CMS top prod 

1.23 

0.55 

0.61 

0.40 

Total 

1.233 

1.264 

1.246 

1.269 


Table 5. Same as Table 4 for the DIS+DY+top fits. 


fits lead to essentially the same value, since in the resummed case the improvement in some 
experiments is compensated by the deterioration of others. 

Turning to the NNLO+NNLL fit results in Table 5, we see that now the effect of 
resummation is more moderate. Effects are small, and in the case again for the fixed-target 
Drell-Yan data, resummation deteriorates the fit quality. Interestingly, the for the LHCb 
Z rapidity data, which, being in the forward region, probe rather large values of x, improves 
substantially with the inclusion of resummation, even at NNLL. Given the small differences 
at the x^ level, we also expect smaller differences at the PDF level, as in the case of the 
DIS-only NNLO+NNLL fit. 

The comparison of the PDFs between the NLO and NLO+NLL DIS+DY+top fits is 
shown in Fig. 10, and the corresponding comparison between the NNLO and NNLO+NNLL 
fits is found in Fig. 11. These can be compared with the corresponding DIS-only fits, see 
Fig. 8 and Fig. 9. In the case of the NLO+NLL fit, the trend is similar to that of the 
DIS-only fit: softer quarks at very large x, and a corresponding enhancement of the valence 
distribution at medium x. At small x, the effect of resummation is negligible as expected. 
We note that both for the total quark singlet and for the total valence PDF the effect of the 
resummation has become more significant than in the DIS-only fit, due to the reduction of 
PDF uncertainties. Indeed, for S(x, for example, the shift of the central value is about 
~ 5% at X ~ 0.5, and the uncertainty bands of the two fits, although still overlapping, are 

NLO+NLL and NNLO+NNLL fits we get ~ 1 Dr the Drell-Yan data. Therefore, the deterioration of 
the x^ of E866 in the resummed fits can be attributed to tension with other datasets, rather than a failure 
of the resummation to correctly describe this dataset. 
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Figure 10. Comparison between the NNPDF3.0 NLO DIS+DY+top fit and the corresponding 
NLO+NLL fit, for as{m%) = 0.118, at a typical LHC scale of = io4 QeV^. 

clearly departing from each other. 

At the NNLO+NNLL level. Fig. 11, the impact of the resummation is as expected even 
smaller. In particular, the PDF uncertainty bands of the fixed order and resummed fits are 
quite similar and they overlap in the entire x range. Large shifts in the central values only 
occur in regions where the PDF uncertainties are large. For example, for the large-x gluon, 
at X ~ 0.3 the resummed central value is ~ 15% smaller than the fixed order one, however 
the PDF uncertainty in this region is substantially larger than the central-value shift. At 
very large x, a similar trend can be seen for the total valence PDF, which however exhibits 
an enhancement at x ~ 0.3, which is as big as the PDF uncertainty. 

4.4 Partonic luminosities 

We now study the impact of the inclusion of threshold resummation on PDFs at the level 
of partonic luminosities. This comparison is useful because it provides direct information 
on how the cross sections for the production of a given final state with invariant mass Mx 
will be affected by the inclusion of resummation in the PDFs. It should be emphasised 
however that in a consistent calculation the impact of the resummation in the PDFs may 
be compensated by a similar sized effect of the resummation in the partonic matrix elements 
of the process under consideration. The consistent comparison of LHC cross sections with 
resummation included both at the PDFs and in the matrix elements is performed in the 
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Figure 11. Same as Fig. 10, now comparing the NNLO DIS+DY+top fit with the corresponding 
NNLO+NNLL fit. 


next section. 

We begin by estimating the effect on the PDF luminosities of the reduced dataset used 
in our baseline hts, as compared to the NNPDF3.0 global ht. Thus in Fig. 12 we compare 
the NNPDF3.0 NLO partonic luminosities for = 0.118, in the global ht and in the 

DIS+DY+top baseline ht. In the upper plots we show the quark-antiquark and quark- 
quark luminosities, and in the lower plots the gluon-gluon and gluon-quark luminosities. 
The calculation has been performed for the LHC 13 TeV, as a function of the mass of the 
hnal state Mx, and results are normalised to the central value of the global ht. 

As we can see in Fig. 12, there are some important differences between the global and 
DIS+DY+top hts. For the qq luminosity, the impact of varying the dataset is small, both 
in terms of central values and of PDF uncertainties, except at very large values of Mx- For 
the qq luminosity, the differences are again only sizeable at large Mx, where the central 
value of the DIS+DY+top ht is softer than that of the global ht, for instance by 10% at 
Mx — 3 TeV. PDF uncertainties are similar in the two cases, and the two hts agree within 
one-sigma. The missing jet data have a stronger impact on the gg and qg luminosities. For 
instance for the gg luminosity above 0.5 TeV, PDF uncertainties increase by a factor two or 
more. Therefore, in order to consistently assess the impact of the resummation, one should 
compare the resummed and hxed-order DIS+DY+top hts, rather than the NNPDF3.0 
global ht, with the resummed hts presented here. 
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LHC 13 TeV, NNPDF3.0 NLO, a3(M^)=0.118 



LHC 13 TeV, NNPDF3.0 NLO, a3(M3)=0.118 
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Figure 12. Comparison of the NNPDF3.0 NLO partonic luminosities for as{rn^) = 0.118, in the 
global fit and in the DIS+DY+top fit which is used as fixed-order baseline for the resummed fits. 
In the upper plots we show the quark-antiquark and quark-quark luminosities, and in the lower 
plots the gluon-gluon and gluon quark luminosities. The calculation has been performed for the 
LHC 13 TeV, as a function of the mass of the final state Mx, and results are normalised to the 
central value of the global fit. 


The comparisons between the DIS+DY+top hxed-order and resummed hts are dis¬ 
played in Figs. 13 (at NLL) and 14 (at NNLL). We see that in all cases the hxed-order 
and resummed hts agree at the level of one sigma, and that the effect of resummation is as 
expected smaller at NNLL than at NLL. In the comparison between NLO and NLO+NLL, 
the qq and qq luminosities are enhanced by about one sigma for Mx ^ 1 TeV, while they 
are suppressed at larger values of Mx- This behaviour follows from the corresponding PDF 
comparisons, where quarks are slightly enhanced at x ~ 0.1 — 0.4 but suppressed for larger 
values of x. This suppression can be sizeable: for Mx — 3 TeV the qq luminosity in the 
NLO+NLL ht is reduced by ~ 15%. In the qq channel, this suppression instead is small 
unless very large values of Mx are probed. The gg and gq luminosities are also suppressed 
at large invariant masses, for instance for gg the suppression is already ~ 10% at 1 TeV, 
though still consistent with the hxed-order ht within the large PDF uncertainties. 

From the corresponding comparison between the NNLO and NNLO+NNLL hts, shown 
in Fig. 14, we see that the effects of resummation are very small everywhere except for the 
largest values of Mx- The central values of the qq, gg and qg luminosities exhibit some 
suppression at very large Mx, but this suppression is not relevant when compared to the 
PDF uncertainties. From the comparison in Fig. 14 we thus conclude that the impact of 
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Figure 13. Same as Fig. 12, now comparing the results of DIS+DY+top fits using either NLO or 
NLO+NLL calculations. 


threshold resummation in a global PDF analysis is only relevant at NLO, while at NNLO it 
appears to be negligible, at least with the current PDF uncertainties. If future data leads 
to substantial reduction of PDF uncertainties at large-x, threshold resummation could be 
relevant even for NNLO fits. 


5 Resummed PDFs: implications for LHC phenomenology 

In this section we discuss the implications of the NLO+NLL and NNLO+NNLL resummed 
PDF sets for LHC phenomenology. Our aim is to quantify, for a variety of processes, the 
difference between using consistently NLO+NLL and NNLO+NNLL calculations at the 
level of both PDFs and matrix elements, and the usual (but inconsistent) approach of using 
resummed partonic cross sections with fixed-order PDFs. 

For illustration, we consider three representative LHC processes for which resummed 
calculations are publicly available, either at the level of total cross sections or of differential 
invariant-mass distributions. We start by considering Higgs production in gluon fusion, 
both for niH = 125 GeV and for a heavy BSM Higgs-like neutral scalar. Note that the 
current recommendation of the Higgs Cross Section Working group for inclusive Higgs 
production in gluon fusion is based on the NNLO+NNLL calculation [35, 173]. We then 
consider threshold resummation for the invariant mass distributions of dileptons in the high- 
mass Drell-Yan process, which is important in many New Physics searches, for example for 
Z' searches. Finally we study the invariant mass distribution of supersymmetric lepton 
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Figure 14. Same as Fig. 13, this time comparing the NNLO and NNLO+NNLL fits. 

(slepton) pair production, a typical final state analysed in electroweak SUSY searches. 
While Higgs production is driven by the gg luminosity, both high-mass Drell-Yan and 
slepton pair production are driven by the quark-antiquark luminosity, which is reasonably 
well constrained even with the reduced dataset used in the present fits. 

A variety of other interesting processes are available in which resummed PDFs should 
be relevant, including top quark differential distributions [62, 63, 174], squark and gluino 
pair production [65, 67] or stop quark pair production [69]. However for most of these 
processes the corresponding resummation codes are not publicly available. 

5.1 SM and BSM Higgs production in gluon fusion 

The accurate calculation of Higgs production via gluon fusion is an essential component 
of the LHC program, since it is required in order to extract Higgs couplings from the 
ATLAS and CMS measurements. As a result of the recent calculation of the inclusive cross 
section at N^LO [24], PDF uncertainties are now one of the dominant theory uncertainties. 
In addition to the characterisation of the SM Higgs boson, many New Physics scenarios 
predict heavy Higgs-like bosons [175-177], and thus it is also important to provide accurate 
predictions for heavy Higgs production for these BSM searches. 

Using the ggHiggs code [121], in Fig. 15 we show the predictions for the total cross 
section for Higgs production in gluon fusion at the LHC 13 TeV, comparing the (N)NLO 
fixed-order results with those of the (N)NLO+(N)NLL resummed calculations, using either 
fixed-order or resummed PDFs. The calculation has been performed in the mtop —t oo limit 
and neglecting finite-width effects, which is sufficient for current purposes. All results are 
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Figure 15. Left: the total cross section for Higgs production in gluon fusion at the LHC 13 TeV 
for different values of the Higgs mass, comparing the predictions of NLO fixed-order with that 
of NLO+NLL resummed calculations (using either fixed order or resummed PDFs). Results are 
shown normalised to the central prediction of the fixed order NLO calculation. Right: the same 
comparison now performed at NNLO. The calculation has been performed using the ggHiggs code. 

normalised to the central value of the fixed-order (N)NLO calculation, and we provide three 
different values of the Higgs mass: mn = 125 GeV, 600 GeV and 2 TeV. 

The comparisons in Fig. 15 are interesting because they show that for the production 
of heavy final states that probe large values of x in the gluon PDF, including resummation 
in the PDFs can cancel out the effect of the resummation in the matrix element. In the 
case of the NLO calculation, the SM Higgs cross section is not affected by resummation of 
the PDFs, but already for mn = 600 GeV, the inclusion of resummed PDFs cancels almost 
half of the enhancement in the hadronic cross section that arises from resummation of the 
matrix element. For an even heavier Higgs, with mn = 2 TeV, the consistent NLO+NLL 
calculation is essentially identical to the NLO result. The trend is similar at NNLO, though 
of course in this case the effect of perturbative corrections beyond the fixed-order NNLO 
calculation is smaller. Note also that PDF uncertainties are substantial at large Higgs 
masses, partly because of the lack of jet data in the baseline and resummed fits. 

Our results demonstrate that using consistently resummed PDFs for SM Higgs pro¬ 
duction at the LHG has no effect, and therefore puts on a more solid ground the current 
HXSWG recommendation, which is based on fixed-order PDFs. This observation is in 
agreement with the findings of Ref. [178] regarding the (lack of) need of N^LO PDFs for 
the SM Higgs production cross section at N^LO. 

5.2 High-Mass Drell-Yan dilepton mass distributions 

At the LHG, high-mass Drell-Yan is one of the most important processes when looking 
for new physics, in particular for new electroweak sectors. For instance, ATLAS and GMS 
have explored a number of BSM signatures in the high-mass tail of neutral-current Drell-Yan 
production [179-181], such as Z' bosons which appear in several new physics scenarios. It 
is therefore interesting to assess the effect of including consistently threshold resummation 
both in the PDFs and in matrix elements, compared to including it only in the matrix 
element while using fixed-order PDFs. 
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Figure 16. Left: dilepton invariant mass distribution for high-mass neutral current Drell-Yan 
production at the LHC 13 TeV, comparing the predictions of fixed-order with that of resummed 
calculations. Results are shown normalised to the central prediction of the fixed-order NLO calcu¬ 
lation. Right: the same comparison at NNLO. 

In Fig. 16 we show the dilepton invariant mass distribution for high-mass neutral current 
Drell-Yan production at the LHC 13 TeV, comparing the predictions of hxed-order and 
resummed calculations. The fixed-order NLO and NNLO predictions have been computed 
with the code Vrap supplemented with threshold resummation as provided by TROLL. In 
Fig. 16 we show the predictions for the dilepton invariant mass distribution, comparing the 
(N)NLO hxed-order results with those of the (N)NLO+(N)NLL resummed calculations, 
using either hxed-order or resummed PDFs. The latter comparison quantihes the mismatch 
when resummed calculations are used with hxed-order PDFs. 

The results are qualitatively consistent with those of the Higgs cross sections in Fig. 15. 
First, we see that even at NLO and at large invariant masses the effect of threshold resum¬ 
mation is moderate: the NLL correction amounts (for hxed-order NLO PDF) to about 
4% at Mu = 2.5 TeV, which is within the current PDF uncertainty. Including the effect 
of the resummation in the PDFs consistently cancels this effect, for example in the range 
for Mu G [1.5, 2.5] TeV the central value of the NLO+NLL calculation agrees with the 
hxed-order NLO result by less than one percent. At NNLO the impact of resummation is 
completely negligible, both at the level of the PDFs and of the matrix elements. 

5.3 Supersymmetric particle production 

The theoretical predictions for high-mass supersymmetric pair production at hadron col¬ 
liders are currently made at NLO, supplemented with either NLL or NNLL resumma¬ 
tion of threshold logarithms. In particular, the NLO+NLL resummed calculations of 
Refs. [182, 183] have been used to produce the benchmark production cross sections at 
^/s = 7 TeV and 13 TeV that are used as by ATLAS and CMS in the theoretical interpre¬ 
tation of their searches for supersymmetry [184-187]. 

An important limitation of these predictions is the mismatch between the hxed-order 
PDFs and the resummed partonic cross sections, which should be more important at high- 
masses, precisely the crucial region for New Physics searches. Thanks to the availability, 
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Figure 17. The NLO+NLL calculation of the invariant mass distribution for slepton pair pro¬ 
duction at the LHC 13 TeV using the Resummino program, using both the NLO and NLO+NLL 
NNPDF3.0 DIS+DY+top PDFs as input. Results are shown as ratios with respect to the NLO 
calculation, using consistently the NLO baseline PDFs. The settings of the SUSY calculation are 
the default ones in Resummino. We use a slepton mass of mj = 564 GeV. 


for the first time, of general-purpose resummed PDFs, it is now possible to consistently 
combine resummed PDFs and matrix elements into a single calculation. It is beyond the 
scope of this work to present a comprehensive study of the impact of NLO+NLL PDFs 
for generic supersymmetric processes. However, for illustrative purposes, in this section we 
will use the public code Resummino [70-72] to compare the effect of resummed PDFs in the 
context of NLO+NLL predictions for electroweak supersymmetric particle pair production 
at the LHC, in particular for slepton pair production. 

Resummino computes resummed and matched predictions for supersymmetric particle 
production at hadron colliders up to the NLO+NLL level. Currently the processes imple¬ 
mented include gaugino-pair production and slepton-pair production. These final states are 
characteristic signatures in electroweak SUSY searches at the LHC [188-191]. Resummino 
is able to compute total cross sections as well as invariant-mass and transverse-momentum 
distributions. In this study we focus on the invariant-mass distribution for slepton pair 
production. Note that the production of sleptons (like many other electroweak SUSY pro¬ 
cesses) is mostly sensitive to the qq luminosity; other processes, such as squark and gluino 
pair production, would be sensitive to other PDF combinations such as qg and gg. 

In Fig. 17 we show the results of the NLO+NLL calculation of the invariant mass 
distribution for slepton pair production at the LHC 13 TeV obtained with Resummino, 
using both the NLO and NLO+NLL NNPDF3.0 DIS+DY+top PDFs as input. Results are 
shown as ratios with respect to the NLO calculation, using consistently the NLO baseline 
PDFs. The settings of the SUSY calculation are the default ones in Resummino. We use a 
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slepton mass of mi = 564 GeV. 

The comparison displayed in Fig. 17 is interesting for a variety of reasons. First of all, we 
see that, nsing the NLO PDFs as inpnt, the NLO+NLL calcnlation (i.e. with resnmmation 
inclnded only in the matrix element) enhances the cross section by several percent, from 2% 
at Mil ~ 1.2 TeV np to 5% at ~ 3 TeV. On the other hand, in the consistent calcnlation 
in which resnmmation is inclnded both at the level of PDFs and of matrix elements, this 
increase is only seen aronnd ~ 1.2 TeV. For higher invariant masses np to Mji ~ 2.5 TeV 
or so the effect of resnmmation in the PDFs cancels the one originating from the matrix 
elements, and the consistent NLO+NLL calcnlation is essentially the same as the NLO 
one. For even higher masses, the NLO+NLL calcnlation is snppressed compared to the 
NLO calcnlation by np to 5% at ~ 3 TeV, thongh in this region PDF nncertainties are 
very large. 

The resnlts of Fig. 17 are consistent with the behavior of the qq Inminosity shown in 
Fig. 13. In particnlar, for an invariant mass of M ~ 1 TeV, the NLO and NLO+NLL PDF 
Inminosities are essentially the same, while for M ~ 3 TeV the NLO+NLL Inminosity is 
snppressed by a factor of approximately 10%, a similar amonnt as that inferred from the 
Resummino plot of Fig. 17. This illnstrates that, for those processes which are dominated by 
a single partonic Inminosity, one can approximately correct a NLO+NLL matrix element 
calcnlation nsing the ratio of PDF Inminosities. This is also consistent with the high-mass 
Drell-Yan resnlts of Fig. 16, which are also driven by the qq Inminosity. 

In snmmary, even thongh we have been able to explore only a limited nnmber of 
resnmmed calcnlations for LHC processes, a consistent trend appears. When the prodnced 
final state has an invariant mass far from threshold, the nse of resnmmed PDFs has a 
rather small effect. However, for heavy final states, the main effect of the resnmmation of 
the PDFs is to compensate the effect of the resnmmation in the matrix element, so that the 
consistent (N)NLO+(N)NLL calcnlation is rather closer to the fixed-order (N)NLO resnlt. 
This shows that nsing resnmmation only in the matrix element bnt not in the PDF can be 
misleading, since it may overestimate cross sections and invariant mass distribntions. This 
is particnlarly the case for NLO+NLL calcnlations, becanse at NNLO+NNLL the effect of 
the resnmmation is mnch smaller, since mnch of it has already been acconnted for in the 
NNLO fixed-order corrections. 

In conclnsion, one shonld in general always nse resnmmed PDFs with threshold re¬ 
snmmed matrix elements. This said, even if the central valne of the consistent (N)NLO+ 
(N)NLL calcnlation is reasonably close to the original fixed-order (N)NLO resnlt, it is in 
general still better to nse the resnmmed calcnlation, since these benefit for instance from 
rednced scale dependence, and thns smaller theoretical nncertainties. 

6 Summary 

In this paper we have presented for the first time global fits of parton distribntions extracted 
at NLO+NLL and NNLO+NNLL accnracy, where the fixed-order partonic cross sections 
have been systematically improved nsing soft-glnon threshold resnmmation. We find that 
the main effect of threshold resnmmation is to snppress the PDFs in the large-x region, as 
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expected given that the fit compensates from the resummation-induced increment in the 
partonic cross sections used in the PDF fit. This suppression is important for all PDF 
flavors for x > 0.1, while at intermediate values of x, 0.01 < x < 0.1, the quark PDFs are 
instead somewhat enhanced due to the sum rules. For smaller values of x, x < 0.01, the 
effect of resummation becomes completely negligible. 

At the level of PDF luminosities at the LHC 13 TeV, we find that at the NLO+NLL 
level the suppression induced by resummation in the PDFs starts to become important for 
Mx ^ 400 GeV in the gg channel, Mx ^ 1 TeV in the qq and qg channels, and Mx ^ 5 TeV 
for the qq channel. The trend is similar at NNLO, but in this case differences between fixed- 
order and resummed PDFs are much smaller. We also find that fixed-order and resummed 
PDFs differ by at most one sigma throughout all the range of Mx- 

We have investigated the corresponding implications at the level of resummed LHC 
cross sections for three different processes: SM and BSM Higgs production in gluon fusion, 
high-mass Drell-Yan pair production and slepton pair production. We find that the effect 
of consistently including resummation in the PDFs can compensate the enhancement from 
resummation of the partonic cross sections, if Mx is large enough. For the production of 
final states with lower Mx, the effect of PDF resummation is negligible. This trend is likely 
to be general: when fitting to data, PDFs adjust to absorb the effect of the resummation in 
the partonic cross sections, and this compensation inevitably persists when extrapolating 
to predictions for new processes. 

Our results illustrate the importance of using the same perturbative order in all the 
components that enter hadronic cross sections: the use of fixed-order PDFs with resummed 
matrix elements can lead to misleading results, especially at high invariant masses, a region 
crucial for new physics searches. The partial cancellation between resummation in PDFs and 
in matrix elements indicates that consistent resummed calculations can be closer to fixed- 
order results. This said, even in the case of a complete cancellation, use of resummation 
would still be advantageous, because of the reduced scale uncertainty. For these reasons, we 
expect that the resummed NNPDF3.0 sets will provide a crucial ingredient, missing so far, 
to improve the precision of all-order resummed calculations for the LHC, and in particular 
those used in searches for new physics. 

The main limitation of the present work is that, since resummed calculations are not 
available for all processes included in the NNPDF3.0 global fit, we have restricted the 
dataset in the resummed fits (and the corresponding baseline) to those processes that can 
be consistently resummed. In particular, we have had to exclude the inclusive jet production 
data and W lepton rapidity distributions. Hence, our resummed sets generate larger PDF 
uncertainties than the NNPDF3.0 global PDFs, especially for gluon-initiated processes. It is 
thus important in the future to provide resummed calculations for these missing processes, 
in order to produce a truly global resummed PDF analysis. 

With the recent start-up of the LHC Run H, the need for precision calculations is even 
more pressing than in Run I, since precision could be the key to uncovering new physics. The 
results of this paper offer for the first time fully consistent threshold resummed calculations, 
which constitute state-of-the-art accuracy for a number of important LHC processes, from 
Higgs to supersymmetric particle production. Therefore, the resummed NNPDF3.0 sets 
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presented here achieve a new milestone in the program of precision phenomenology at the 
LHC. 


Delivery 

The resummed calculations for DIS structure function and Drell-Yan distributions used in 
this work have been obtained with the new code TROLL, version v3.0. This code is publicly 
available from 


http: //www.ge. infn. it/~bonvini/troll/ 

and can be used by any interested parties to compute their own predictions for resummed 
observables. 

Our resummed sets are available in the LHAPDF6 format from the authors upon request. 
The available sets are: 


NNPDF30_nlo_disdytop 

NNPDF30_nnlo_disdytop 

NNPDF30_nll_disdytop 

NNPDF30_nnll_disdytop 

which stand for the NLO and NNLO baseline fits, and their NLO+NLL and NNLO+NNLL 
resummed counterparts. All these PDF sets are provided for = 0.118 with a 

maximum of ny = 5 active flavors. 
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